Relation between HLA and copy number variation of steroid 21-hydroxylase in a Swedish cohort of patients with autoimmune Addison's disease.

OBJECTIVE: Autoantibodies against the adrenal enzyme 21-hydroxylase is a hallmark manifestation in autoimmune Addison's disease (AAD). Steroid 21-hydroxylase is encoded by CYP21A2, which is located in the human leucocyte antigen (HLA) region together with the highly similar pseudogene CYP21A1P. A high level of copy number variation is seen for the 2 genes, and therefore, we asked whether genetic variation of the CYP21 genes is associated with AAD. DESIGN: Case-control study on patients with AAD and healthy controls. METHODS: Using next-generation DNA sequencing, we estimated the copy number of CYP21A2 and CYP21A1P, together with HLA alleles, in 479 Swedish patients with AAD and autoantibodies against 21-hydroxylase and in 1393 healthy controls. RESULTS: With 95% of individuals carrying 2 functional 21-hydroxylase genes, no difference in CYP21A2 copy number was found when comparing patients and controls. In contrast, we discovered a lower copy number of the pseudogene CYP21A1P among AAD patients (P = 5 × 10-44), together with associations of additional nucleotide variants, in the CYP21 region. However, the strongest association was found for HLA-DQB1*02:01 (P = 9 × 10-63), which, in combination with the DRB1*04:04-DQB1*03:02 haplotype, imposed the greatest risk of AAD. CONCLUSIONS: We identified strong associations between copy number variants in the CYP21 region and risk of AAD, although these associations most likely are due to linkage disequilibrium with disease-associated HLA class II alleles.

Identification and functional characterization of a novel susceptibility locus for small vessel vasculitis with MPO-ANCA.

OBJECTIVE: To identify and characterize genetic loci associated with the risk of developing ANCA-associated vasculitides (AAV). METHODS: Genetic association analyses were performed after Illumina sequencing of 1853 genes and subsequent replication with genotyping of selected single nucleotide polymorphisms in a total cohort of 1110 Scandinavian cases with granulomatosis with polyangiitis or microscopic polyangiitis, and 1589 controls. A novel AAV-associated single nucleotide polymorphism was analysed for allele-specific effects on gene expression using luciferase reporter assay. RESULTS: PR3-ANCA+ AAV was significantly associated with two independent loci in the HLA-DPB1/HLA-DPA1 region [rs1042335, P = 6.3 × 10-61, odds ratio (OR) 0.10; rs9277341, P = 1.5 × 10-44, OR 0.22] and with rs28929474 in the SERPINA1 gene (P = 2.7 × 10-10, OR 2.9). MPO-ANCA+ AAV was significantly associated with the HLA-DQB1/HLA-DQA2 locus (rs9274619, P = 5.4 × 10-25, OR 3.7) and with a rare variant in the BACH2 gene (rs78275221, P = 7.9 × 10-7, OR 3.0), the latter a novel susceptibility locus for MPO-ANCA+ granulomatosis with polyangiitis/microscopic polyangiitis. The rs78275221-A risk allele reduced luciferase gene expression in endothelial cells, specifically, as compared with the non-risk allele. CONCLUSION: We identified a novel susceptibility locus for MPO-ANCA+ AAV and propose that the associated variant is of mechanistic importance, exerting a regulatory function on gene expression in specific cell types.

A novel canine reference genome resolves genomic architecture and uncovers transcript complexity.

We present GSD_1.0, a high-quality domestic dog reference genome with chromosome length scaffolds and contiguity increased 55-fold over CanFam3.1. Annotation with generated and existing long and short read RNA-seq, miRNA-seq and ATAC-seq, revealed that 32.1% of lifted over CanFam3.1 gaps harboured previously hidden functional elements, including promoters, genes and miRNAs in GSD_1.0. A catalogue of canine "dark" regions was made to facilitate mapping rescue. Alignment in these regions is difficult, but we demonstrate that they harbour trait-associated variation. Key genomic regions were completed, including the Dog Leucocyte Antigen (DLA), T Cell Receptor (TCR) and 366 COSMIC cancer genes. 10x linked-read sequencing of 27 dogs (19 breeds) uncovered 22.1 million SNPs, indels and larger structural variants. Subsequent intersection with protein coding genes showed that 1.4% of these could directly influence gene products, and so provide a source of normal or aberrant phenotypic modifications.

Whole-genome genotyping and resequencing reveal the association of a deletion in the complex interferon alpha gene cluster with hypothyroidism in dogs.

BACKGROUND: Hypothyroidism is a common complex endocrinopathy that typically has an autoimmune etiology, and it affects both humans and dogs. Genetic and environmental factors are both known to play important roles in the disease development. In this study, we sought to identify the genetic risk factors potentially involved in the susceptibility to the disease in the high-risk Giant Schnauzer dog breed. RESULTS: By employing genome-wide association followed by fine-mapping (top variant p-value = 5.7 × 10- 6), integrated with whole-genome resequencing and copy number variation analysis, we detected a ~ 8.9 kbp deletion strongly associated (p-value = 0.0001) with protection against development of hypothyroidism. The deletion is located between two predicted Interferon alpha (IFNA) genes and it may eliminate functional elements potentially involved in the transcriptional regulation of these genes. Remarkably, type I IFNs have been extensively associated to human autoimmune hypothyroidism and general autoimmunity. Nonetheless, the extreme genomic complexity of the associated region on CFA11 warrants further long-read sequencing and annotation efforts in order to ascribe functions to the identified deletion and to characterize the canine IFNA gene cluster in more detail. CONCLUSIONS: Our results expand the current knowledge on genetic determinants of canine hypothyroidism by revealing a significant link with the human counterpart disease, potentially translating into better diagnostic tools across species, and may contribute to improved canine breeding strategies.

A rare regulatory variant in the MEF2D gene affects gene regulation and splicing and is associated with a SLE sub-phenotype in Swedish cohorts.

Systemic lupus erythematosus (SLE) is an autoimmune disorder with heterogeneous clinical presentation and complex etiology involving the interplay between genetic, epigenetic, environmental and hormonal factors. Many common SNPs identified by genome wide-association studies (GWAS) explain only a small part of the disease heritability suggesting the contribution from rare genetic variants, undetectable in GWAS, and complex epistatic interactions. Using targeted re-sequencing of coding and conserved regulatory regions within and around 215 candidate genes selected on the basis of their known role in autoimmunity and genes associated with canine immune-mediated diseases, we identified a rare regulatory variant rs200395694:G > T located in intron 4 of the MEF2D gene encoding the myocyte-specific enhancer factor 2D transcription factor and associated with SLE in Swedish cohorts (504 SLE patients and 839 healthy controls, p = 0.014, CI = 1.1-10). Fisher's exact test revealed an association between the genetic variant and a triad of disease manifestations including Raynaud, anti-U1-ribonucleoprotein (anti-RNP), and anti-Smith (anti-Sm) antibodies (p = 0.00037) among the patients. The DNA-binding activity of the allele was further studied by EMSA, reporter assays, and minigenes. The region has properties of an active cell-specific enhancer, differentially affected by the alleles of rs200395694:G > T. In addition, the risk allele exerts an inhibitory effect on the splicing of the alternative tissue-specific isoform, and thus may modify the target gene set regulated by this isoform. These findings emphasize the potential of dissecting traits of complex diseases and correlating them with rare risk alleles with strong biological effects.

Common genetic variation in the autoimmune regulator (AIRE) locus is associated with autoimmune Addison's disease in Sweden.

Autoimmune Addison's disease (AAD) is the predominating cause of primary adrenal failure. Despite its high heritability, the rarity of disease has long made candidate-gene studies the only feasible methodology for genetic studies. Here we conducted a comprehensive reinvestigation of suggested AAD risk loci and more than 1800 candidate genes with associated regulatory elements in 479 patients with AAD and 2394 controls. Our analysis enabled us to replicate many risk variants, but several other previously suggested risk variants failed confirmation. By exploring the full set of 1800 candidate genes, we further identified common variation in the autoimmune regulator (AIRE) as a novel risk locus associated to sporadic AAD in our study. Our findings not only confirm that multiple loci are associated with disease risk, but also show to what extent the multiple risk loci jointly associate to AAD. In total, risk loci discovered to date only explain about 7% of variance in liability to AAD in our study population.

Cytokine Autoantibody Screening in the Swedish Addison Registry Identifies Patients With Undiagnosed APS1.

Context: Autoimmune polyendocrine syndrome type 1 (APS1) is a monogenic disorder that features autoimmune Addison disease as a major component. Although APS1 accounts for only a small fraction of all patients with Addison disease, early identification of these individuals is vital to prevent the potentially lethal complications of APS1. Objective: To determine whether available serological and genetic markers are valuable screening tools for the identification of APS1 among patients diagnosed with Addison disease. Design: We systematically screened 677 patients with Addison disease enrolled in the Swedish Addison Registry for autoantibodies against interleukin-22 and interferon-α4. Autoantibody-positive patients were investigated for clinical manifestations of APS1, additional APS1-specific autoantibodies, and DNA sequence and copy number variations of AIRE. Results: In total, 17 patients (2.5%) displayed autoantibodies against interleukin-22 and/or interferon-α4, of which nine were known APS1 cases. Four patients previously undiagnosed with APS1 fulfilled clinical, genetic, and serological criteria. Hence, we identified four patients with undiagnosed APS1 with this screening procedure. Conclusion: We propose that patients with Addison disease should be routinely screened for cytokine autoantibodies. Clinical or serological support for APS1 should warrant DNA sequencing and copy number analysis of AIRE to enable early diagnosis and prevention of lethal complications.

Multiple regulatory variants located in cell type-specific enhancers within the PKP2 locus form major risk and protective haplotypes for canine atopic dermatitis in German shepherd dogs.

BACKGROUND: Canine atopic dermatitis (CAD) is a chronic inflammatory skin disease triggered by allergic reactions involving IgE antibodies directed towards environmental allergens. We previously identified a ~1.5 Mb locus on canine chromosome 27 associated with CAD in German shepherd dogs (GSDs). Fine-mapping indicated association closest to the PKP2 gene encoding plakophilin 2. RESULTS: Additional genotyping and association analyses in GSDs combined with control dogs from five breeds with low-risk for CAD revealed the top SNP 27:19,086,778 (p = 1.4 × 10(-7)) and a rare ~48 kb risk haplotype overlapping the PKP2 gene and shared only with other high-risk CAD breeds. We selected altogether nine SNPs (four top-associated in GSDs and five within the ~48 kb risk haplotype) that spanned ~280 kb forming one risk haplotype carried by 35 % of the GSD cases and 10 % of the GSD controls (OR = 5.1, p = 5.9 × 10(-5)), and another haplotype present in 85 % of the GSD cases and 98 % of the GSD controls and conferring a protective effect against CAD in GSDs (OR = 0.14, p = 0.0032). Eight of these SNPs were analyzed for transcriptional regulation using reporter assays where all tested regions exerted regulatory effects on transcription in epithelial and/or immune cell lines, and seven SNPs showed allelic differences. The DNA fragment with the top-associated SNP 27:19,086,778 displayed the highest activity in keratinocytes with 11-fold induction of transcription by the risk allele versus 8-fold by the control allele (pdifference = 0.003), and also mapped close (~3 kb) to an ENCODE skin-specific enhancer region. CONCLUSIONS: Our experiments indicate that multiple CAD-associated genetic variants located in cell type-specific enhancers are involved in gene regulation in different cells and tissues. No single causative variant alone, but rather multiple variants combined in a risk haplotype likely contribute to an altered expression of the PKP2 gene, and possibly nearby genes, in immune and epithelial cells, and predispose GSDs to CAD.

Utilizing the Dog Genome in the Search for Novel Candidate Genes Involved in Glioma Development-Genome Wide Association Mapping followed by Targeted Massive Parallel Sequencing Identifies a Strongly Associated Locus.

Gliomas are the most common form of malignant primary brain tumors in humans and second most common in dogs, occurring with similar frequencies in both species. Dogs are valuable spontaneous models of human complex diseases including cancers and may provide insight into disease susceptibility and oncogenesis. Several brachycephalic breeds such as Boxer, Bulldog and Boston Terrier have an elevated risk of developing glioma, but others, including Pug and Pekingese, are not at higher risk. To identify glioma-associated genetic susceptibility factors, an across-breed genome-wide association study (GWAS) was performed on 39 dog glioma cases and 141 controls from 25 dog breeds, identifying a genome-wide significant locus on canine chromosome (CFA) 26 (p = 2.8 x 10-8). Targeted re-sequencing of the 3.4 Mb candidate region was performed, followed by genotyping of the 56 SNVs that best fit the association pattern between the re-sequenced cases and controls. We identified three candidate genes that were highly associated with glioma susceptibility: CAMKK2, P2RX7 and DENR. CAMKK2 showed reduced expression in both canine and human brain tumors, and a non-synonymous variant in P2RX7, previously demonstrated to have a 50% decrease in receptor function, was also associated with disease. Thus, one or more of these genes appear to affect glioma susceptibility.

A Multi-Breed Genome-Wide Association Analysis for Canine Hypothyroidism Identifies a Shared Major Risk Locus on CFA12.

Hypothyroidism is a complex clinical condition found in both humans and dogs, thought to be caused by a combination of genetic and environmental factors. In this study we present a multi-breed analysis of predisposing genetic risk factors for hypothyroidism in dogs using three high-risk breeds--the Gordon Setter, Hovawart and the Rhodesian Ridgeback. Using a genome-wide association approach and meta-analysis, we identified a major hypothyroidism risk locus shared by these breeds on chromosome 12 (p = 2.1x10(-11)). Further characterisation of the candidate region revealed a shared ~167 kb risk haplotype (4,915,018-5,081,823 bp), tagged by two SNPs in almost complete linkage disequilibrium. This breed-shared risk haplotype includes three genes (LHFPL5, SRPK1 and SLC26A8) and does not extend to the dog leukocyte antigen (DLA) class II gene cluster located in the vicinity. These three genes have not been identified as candidate genes for hypothyroid disease previously, but have functions that could potentially contribute to the development of the disease. Our results implicate the potential involvement of novel genes and pathways for the development of canine hypothyroidism, raising new possibilities for screening, breeding programmes and treatments in dogs. This study may also contribute to our understanding of the genetic etiology of human hypothyroid disease, which is one of the most common endocrine disorders in humans.

Multiple Changes of Gene Expression and Function Reveal Genomic and Phenotypic Complexity in SLE-like Disease.

The complexity of clinical manifestations commonly observed in autoimmune disorders poses a major challenge to genetic studies of such diseases. Systemic lupus erythematosus (SLE) affects humans as well as other mammals, and is characterized by the presence of antinuclear antibodies (ANA) in patients' sera and multiple disparate clinical features. Here we present evidence that particular sub-phenotypes of canine SLE-related disease, based on homogenous (ANA(H)) and speckled ANA (ANA(S)) staining pattern, and also steroid-responsive meningitis-arteritis (SRMA) are associated with different but overlapping sets of genes. In addition to association to certain MHC alleles and haplotypes, we identified 11 genes (WFDC3, HOMER2, VRK1, PTPN3, WHAMM, BANK1, AP3B2, DAPP1, LAMTOR3, DDIT4L and PPP3CA) located on five chromosomes that contain multiple risk haplotypes correlated with gene expression and disease sub-phenotypes in an intricate manner. Intriguingly, the association of BANK1 with both human and canine SLE appears to lead to similar changes in gene expression levels in both species. Our results suggest that molecular definition may help unravel the mechanisms of different clinical features common between and specific to various autoimmune disease phenotypes in dogs and humans.

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds.

BACKGROUND: The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants. RESULTS: We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves. CONCLUSIONS: We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Constitutive activation of the ERK pathway in melanoma and skin melanocytes in Grey horses.

BACKGROUND: Constitutive activation of the ERK pathway, occurring in the vast majority of melanocytic neoplasms, has a pivotal role in melanoma development. Different mechanisms underlie this activation in different tumour settings. The Grey phenotype in horses, caused by a 4.6 kb duplication in intron 6 of Syntaxin 17 (STX17), is associated with a very high incidence of cutaneous melanoma, but the molecular mechanism behind the melanomagenesis remains unknown. Here, we investigated the involvement of the ERK pathway in melanoma development in Grey horses. METHODS: Grey horse melanoma tumours, cell lines and normal skin melanocytes were analyzed with help of indirect immunofluorescence and immunoblotting for the expression of phospho-ERK1/2 in comparison to that in non-grey horse and human counterparts. The mutational status of BRAF, RAS, GNAQ, GNA11 and KIT genes in Grey horse melanomas was determined by direct sequencing. The effect of RAS, RAF and PI3K/AKT pathways on the activation of the ERK signaling in Grey horse melanoma cells was investigated with help of specific inhibitors and immunoblotting. Individual roles of RAF and RAS kinases on the ERK activation were examined using si-RNA based approach and immunoblotting. RESULTS: We found that the ERK pathway is constitutively activated in Grey horse melanoma tumours and cell lines in the absence of somatic activating mutations in BRAF, RAS, GNAQ, GNA11 and KIT genes or alterations in the expression of the main components of the pathway. The pathway is mitogenic and is mediated by BRAF, CRAF and KRAS kinases. Importantly, we found high activation of the ERK pathway also in epidermal melanocytes, suggesting a general predisposition to melanomagenesis in these horses. CONCLUSIONS: These findings demonstrate that the presence of the intronic 4.6 kb duplication in STX17 is strongly associated with constitutive activation of the ERK pathway in melanocytic cells in Grey horses in the absence of somatic mutations commonly linked to the activation of this pathway during melanomagenesis. These findings are consistent with the universal importance of the ERK pathway in melanomagenesis and may have valuable implications for human melanoma research.

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication.

The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

A simple repeat polymorphism in the MITF-M promoter is a key regulator of white spotting in dogs.

The white spotting locus (S) in dogs is colocalized with the MITF (microphtalmia-associated transcription factor) gene. The phenotypic effects of the four S alleles range from solid colour (S) to extreme white spotting (s(w)). We have investigated four candidate mutations associated with the s(w) allele, a SINE insertion, a SNP at a conserved site and a simple repeat polymorphism all associated with the MITF-M promoter as well as a 12 base pair deletion in exon 1B. The variants associated with white spotting at all four loci were also found among wolves and we conclude that none of these could be a sole causal mutation, at least not for extreme white spotting. We propose that the three canine white spotting alleles are not caused by three independent mutations but represent haplotype effects due to different combinations of causal polymorphisms. The simple repeat polymorphism showed extensive diversity both in dogs and wolves, and allele-sharing was common between wolves and white spotted dogs but was non-existent between solid and spotted dogs as well as between wolves and solid dogs. This finding was unexpected as Solid is assumed to be the wild-type allele. The data indicate that the simple repeat polymorphism has been a target for selection during dog domestication and breed formation. We also evaluated the significance of the three MITF-M associated polymorphisms with a Luciferase assay, and found conclusive evidence that the simple repeat polymorphism affects promoter activity. Three alleles associated with white spotting gave consistently lower promoter activity compared with the allele associated with solid colour. We propose that the simple repeat polymorphism affects cooperativity between transcription factors binding on either flanking sides of the repeat. Thus, both genetic and functional evidence show that the simple repeat polymorphism is a key regulator of white spotting in dogs.

Thorough investigation of a canine autoinflammatory disease (AID) confirms one main risk locus and suggests a modifier locus for amyloidosis.

Autoinflammatory disease (AID) manifests from the dysregulation of the innate immune system and is characterised by systemic and persistent inflammation. Clinical heterogeneity leads to patients presenting with one or a spectrum of phenotypic signs, leading to difficult diagnoses in the absence of a clear genetic cause. We used separate genome-wide SNP analyses to investigate five signs of AID (recurrent fever, arthritis, breed specific secondary dermatitis, otitis and systemic reactive amyloidosis) in a canine comparative model, the pure bred Chinese Shar-Pei. Analysis of 255 DNA samples revealed a shared locus on chromosome 13 spanning two peaks of association. A three-marker haplotype based on the most significant SNP (p<2.6×10(-8)) from each analysis showed that one haplotypic pair (H13-11) was present in the majority of AID individuals, implicating this as a shared risk factor for all phenotypes. We also noted that a genetic signature (F ST) distinguishing the phenotypic extremes of the breed specific Chinese Shar-Pei thick and wrinkled skin, flanked the chromosome 13 AID locus; suggesting that breed development and differentiation has played a parallel role in the genetics of breed fitness. Intriguingly, a potential modifier locus for amyloidosis was revealed on chromosome 14, and an investigation of candidate genes from both this and the chromosome 13 regions revealed significant (p<0.05) renal differential expression in four genes previously implicated in kidney or immune health (AOAH, ELMO1, HAS2 and IL6). These results illustrate that phenotypic heterogeneity need not be a reflection of genetic heterogeneity, and that genetic modifiers of disease could be masked if syndromes were not first considered as individual clinical signs and then as a sum of their component parts.

Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis.

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

Complex inheritance of melanoma and pigmentation of coat and skin in Grey horses.

The dominant phenotype of greying with age in horses, caused by a 4.6-kb duplication in intron 6 of STX17, is associated with a high incidence of melanoma and vitiligo-like skin depigmentation. However, the progressive greying and the incidence of melanoma, vitiligo-like depigmentation, and amount of speckling in these horses do not follow a simple inheritance pattern. To understand their inheritance, we analysed the melanoma grade, grey level, vitiligo grade, and speckling grade of 1,119 Grey horses (7,146 measurements) measured in six countries over a 9-year period. We estimated narrow sense heritability (h(2)), and we decomposed this parameter into polygenic heritability (h(2) (POLY)), heritability due to the Grey (STX17) mutation (h(2) (STX17)), and heritability due to agouti (ASIP) locus (h(2) (ASIP)). A high heritability was found for greying (h(2) = 0.79), vitiligo (h(2) = 0.63), and speckling (h(2) = 0.66), while a moderate heritability was estimated for melanoma (h(2) = 0.37). The additive component of ASIP was significantly different from zero only for melanoma (h(2) (ASIP) = 0.02). STX17 controlled large proportions of phenotypic variance (h(2) (STX17) = 0.18-0.55) and overall heritability (h(2) (STX17)/h(2) = 0.28-0.83) for all traits. Genetic correlations among traits were estimated as moderate to high, primarily due to the effects of the STX17 locus. Nevertheless, the correlation between progressive greying and vitiligo-like depigmentation remained large even after taking into account the effects of STX17. We presented a model where four traits with complex inheritance patterns are strongly influenced by a single mutation. This is in line with evidence of recent studies in domestic animals indicating that some complex traits are, in addition to the large number of genes with small additive effects, influenced by genes of moderate-to-large effect. Furthermore, we demonstrated that the STX17 mutation explains to a large extent the moderate to high genetic correlations among traits, providing an example of strong pleiotropic effects caused by a single gene.

Copy number expansion of the STX17 duplication in melanoma tissue from Grey horses.

BACKGROUND: Greying with age in horses is an autosomal dominant trait, associated with loss of hair pigmentation, melanoma and vitiligo-like depigmentation. We recently identified a 4.6 kb duplication in STX17 to be associated with the phenotype. The aims of this study were to investigate if the duplication in Grey horses shows copy number variation and to exclude that any other polymorphism is uniquely associated with the Grey mutation. RESULTS: We found little evidence for copy number expansion of the duplicated sequence in blood DNA from Grey horses. In contrast, clear evidence for copy number expansions was indicated in five out of eight tested melanoma tissues or melanoma cell lines. A tendency of a higher copy number in aggressive tumours was also found. Massively parallel resequencing of the ~350 kb Grey haplotype did not reveal any additional mutations perfectly associated with the phenotype, confirming the duplication as the true causative mutation. We identified three SNP alleles that were present in a subset of Grey haplotypes within the 350 kb region that shows complete linkage disequilibrium with the causative mutation. Thus, these three nucleotide substitutions must have occurred subsequent to the duplication, consistent with our interpretation that the Grey mutation arose more than 2,000 years before present. CONCLUSIONS: These results suggest that the mutation acts as a melanoma-driving regulatory element. The elucidation of the mechanistic features of the duplication will be of considerable interest for the characterization of these horse melanomas as well as for the field of human melanoma research.